arXiv:1506.01405vl [astro-ph.CO] 3 Jun 2015 


Mon. Not. R. Astron. Soc. 000, 1-11 (2015) 


Printed 5 June 2015 (MN style file v2.2) 


Major Mergers Going Notts: Challenges for Modern Halo Finders 

Peter Behroozi\* Alexander Knebe^’^, Frazer R. Pearce^, Pascal Elahi^, Jiaxin Han^, 
Hanni Lux^’^, Yao-Yuan Mao^’^, Stuart 1. Muldrew^^, Doug Potter^Chaichalit Srisawat^^ 

^ Space Telescope Science Institute, Baltimore, MD 21218, USA 

^Departamento de Fi'sica Teorica, Modulo 15, Facultad de Ciencias, Universidad Autonoma de Madrid (UAM), 28049 Madrid, Spain 
^Astro-UAM, UAM, Unidad Asociada CSIC 

‘^School of Physics & Astronomy, University of Nottingham, Nottingham NG7 2RD, UK 
^ Sydney Institute for Astronomy, School of Physics, A28, The University of Sydney, NSW 2006, Australia 
^Institute for Computational Cosmology, Department of Physics, Durham University, South Road, Durham DHl 3LE, UK 
Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford, 0X1 3RH, UK 
^Kavli Institute for Particle Astrophysics and Cosmology & Physics Department, Stanford University, Stanford, CA 94305, USA 
'^SLAC National Accelerator Laboratory, Menlo Park, CA, 94025, USA 
Department of Physics and Astronomy, University of Leicester, University Road, Leicester, LEI 7R11, UK 
Institute for Computational Science, University of Zurich, 8057 Zurich, Switzerland 
Department of Physics and Astronomy, University of Sussex, Brighton BNl 9QH, UK 


Released 5 June 2015 


ABSTRACT 

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo 
finders. We compare five halo finding algorithms’ (AHF, HBT, Rockstar, SubFind, and 
VELOCIraptor) recovery of halo properties for both isolated and cosmological major 
mergers. We find that halo positions and velocities are often robust, but mass biases exist for 
every technique. The algorithms also show strong disagreement in the prevalence and dura¬ 
tion of major mergers, especially at high redshifts (z > 1). This raises significant uncertainties 
for theoretical models that require major mergers for, e.g., galaxy morphology changes, size 
changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not 
using temporal information also show host halo and subhalo relationship swaps over succes¬ 
sive timesteps, requiring careful merger tree construction to avoid problematic mass accretion 
histories. We suggest that future algorithms should combine phase-space and temporal infor¬ 
mation to avoid the issues presented. 

Key words: methods:numerical - dark matter - galaxies: haloes 


1 INTRODUCTION 

In the Lambda Cold Dark Matter paradigm, nonlinear gravitational 
collapse of matter overdensities yields self-bound structures known 
as “haloes.” Smaller haloes merge onto larger ones continuously, 
and are called “subhaloes” as long as they remain distinguishable 
within the radius of the larger halo. Major mergers - i.e., merg¬ 
ers between two haloes of similar mass - occur rarely. Estimates 
from recent simulations suggest that haloes at z = 0 experience 
(on average) one merger per halo per 10 Gyrs (Fakhouri & Ma 
2008; Wetzel, Cohn & White 2009; Behroozi, Wechsler & Conroy 
2013b). 

Despite their infrequency, major mergers have been invoked 
to explain a surprisingly broad range of galaxy phenomena. Galaxy 
growth correlates tightly with halo growth (see Leauthaud et al. 
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2012b; Wang et al. 2013; Behroozi, Wechsler & Conroy 2013a,b; 
Moster, Naab & White 2013; Behroozi & Silk 2015, for recent 
constraints), so the significant disturbance to haloes in major 
mergers could also imply significant changes in observable galaxy 
properties. Merger-linked phenomena with significant recent inter¬ 
est include active galactic nuclei (AGN) activity (Kocevski et al. 
2012; Newton & Kay 2013) and associated black hole growth 
(Treister et al. 2012; Bonoli, Mayer & Callegari 2014), Ultra- 
Luminous InfraRed Galaxy (ULIRG) triggering (Kartaltepe et al. 
2010; Draper & Ballantyne 2012), galaxy morphology and size 
changes (Bernardi et al. 2011; Prieto et al. 2013), galaxy number 
density changes (Lotz et al. 2011; Behroozi et al. 2013a), velocity 
dispersion evolution (Oser et al. 2012), star formation quench¬ 
ing/triggering (Kaviraj et al. 2013; Behroozi et al. 2015), galactic 
winds (Hopkins et al. 2013; Rupke & Veilleux 2013), buildup of 
intracluster light (Laporte et al. 2013), buildup of spheroidal bulges 
(Sales et al. 2012; Wilman et al. 2013), dispersal of magnetic fields 
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(Xu et al. 2010), and creation of tidal shells (Wang et al. 2012). 
They also represent an important systematic for cluster analysis, 
including violations of hydrostatic equilibrium for X-ray masses 
(Akahori & Yoshikawa 2010; Takizawa, Nagino & Matsushita 
2010; Bourdin et al. 2011; Nelson et al. 2012), biases in Sunyaev- 
ZeTdovich signals (Rudd & Nagai 2009; AMI Consortium et al. 
2011; Krause et al. 2012), and incidence of Bullet Cluster-like 
systems (Thompson, Dave & Nagamine 2014; Bouillot et al. 
2015). 

Predicting how major mergers will impact observables of¬ 
ten involves a dark matter or hydrodynamical simulation (see 
Kuhlen, Vogelsberger & Angulo 2012, for a review), a halo finder 
to convert the simulation particle data into a list of haloes and 
their properties (see Knebe et al. 2011, 2013b, for reviews), a 
merger tree algorithm to connect haloes across redshifts (see 
Srisawat et al. 2013, for a review), and optionally a theoreti¬ 
cal model for galaxy formation (see Somerville & Dave 2014, 
for a review). The role of the halo finder in this process has 
been investigated in a recent series of workshops (including 
“Haloes Going MAD,” “Subhaloes Going Notts,” and “Sussing 
Merger Trees”) and papers (Knebe et al. 2011; Onions et al. 2012; 
Elahi et al. 2013; Onions et al. 2013; Pujol et al. 2014; Knebe et al. 
2013a,b; Srisawat et al. 2013; Avila et al. 2014; Lee et al. 2014; 
Hoffmann et al. 2014). This paper, arising out of the “Subhaloes 
Going Notts” and “Sussing Merger Trees” workshops, continues 
this pattern with an investigation into how halo finders treat major 
mergers. 

Halo finder recovery of very minor subhaloes (mass ratios 
<1:10) has already been investigated (Muldrew, Pearce & Power 
2011; Knebe et al. 2011, 2013b; Onions et al. 2012). Finders which 
use particle positions alone to initially classify subhaloes are able 
to perform just as well as finders which use additional information 
(e.g., particle velocities or historical positions) in the outer halo, as 
long as gravitational unbinding is performed (Onions et al. 2012; 
Knebe et al. 2013b). As the larger “host” halo has a much larger 
velocity dispersion than the smaller subhalo, particle binding en¬ 
ergies can very effectively distinguish between particles belonging 
to the host halo and to the subhalo. In major mergers (mass ratios 
greater than 1:3), the host and the subhalo have similar velocity 
dispersions, making particle assignment much harder. Additionally, 
the choice of which halo to call the “host halo” and which to call 
the “subhalo” can be ambiguous for major mergers, and can change 
over time unless temporal information is used (Tweed et al. 2009; 
Han et al. 2012; Behroozi, Wechsler & Wu 2013; Srisawat et al. 
2013). Hence, we investigate halo finders’ abilities not only to re¬ 
cover halo properties in major mergers but also to follow halo prop¬ 
erties smoothly across simulation timesteps. 

We divide the results into several sections. In §2, we briefly 
describe the participating halo finders. We describe “static” tests 
of halo finding, with overlapping mock Navarro, Frenk & White 
(1997) (NFW) profiles, in §3. “Dynamical” tests are presented in 
§4, where two mock NFW profiles are allowed to merge in an iso¬ 
lated simulation. We present tests drawn from cosmological sim¬ 
ulations in §5. Finally, we discuss the impact of these results in 
§6 and summarize our conclusions in §7. Throughout this work, 
halo masses are calculated as spherical overdensities, and host halo 
masses include all substructure masses. 


2 COMMON TERMS AND HALO FINDER 
DESCRIPTIONS 

In this section, we define common terms and briefly describe the 
participating halo finders (AHF, HBT, Rockstar, SubFind, and 
VELOCIraptor). These descriptions include the overall algo¬ 
rithm employed and particle information used (e.g., positions, ve¬ 
locities, halo membership at previous snapshots), the method for 
assigning particles in major mergers, and the method for deciding 
which of two overlapping haloes is the host halo in major mergers. 
Names following the halo finders are the co-authors who ran the 
halo finders for this study and provided the following descriptions, 
who are not always the same as the original halo finder authors. 


2.1 Common Terms 

Throughout this paper, pc = hlP j {SnG) refers to the critical den¬ 
sity. We use Ryyyc to indicate the radius from a halo centre within 
which the average enclosed density is YYY x Pc (including sub¬ 
structure). Similarly, Myyyc refers to the total mass enclosed within 
Ryyyc- Some halo finders also use Rvii^ corresponding to an aver¬ 
age enclosed density as defined in Bryan & Norman (1998). 
The term “Vmax” refers to the maximum circular velocity; i.e., the 
maximum value of GM{< R)/R over a halo’s radial mass profile. 
Finally, “position-space” information refers to particle positions, 
“velocity-space” information refers to particle velocities, “phase- 
space” information refers to both particle positions and velocities, 
and “temporal” information refers to the evolution of particles’ halo 
memberships over time. 


2.2 AHF (Knebe) 

The halo finder AHF* (AMIGA Halo Finder, Knollmann & Knebe 
2009), is an improvement of the MHF halo finder 
(Gill, Knebe & Gibson 2004), which employs a recursively 
refined grid to locate local overdensities in the density field. The 
identified density peaks are then treated as centres of prospective 
haloes. The resulting grid hierarchy is further utilized to generate 
a halo tree readily containing the information which halo is a 
(prospective) host and subhalo, respectively. Halo properties are 
calculated based on the list of particles asserted to be gravita¬ 
tionally bound to the respective density peak. To generate this 
list of particles we employ an iterative procedure starting from 
an initial guess of particles. This initial guess is based upon the 
distance of each prospective centre to its nearest more massive 
(sub-)halo where all particles within a sphere of radius half this 
distance are considered prospective (sub-)halo constituents. This 
tentative particle list is then used in an iterative procedure to 
remove unbound particles and the final particle list is truncated at 
some user pre-defined overdensity criterion. 

The tree for each halo consists of one trunk and several 
branches where the trunk is the continuation of the main host 
halo and the branches represent the subhaloes (see Fig. 1 in 
Knollmann & Knebe 2009). While there are various options in 
AHF to pick the trunk, the default mode (also applied here) is to 
recursively follow the branch containing the most particles. This 
choice certainly leaves its imprint during major merger events stud¬ 
ied here. 


* AHF is freely available from http : / /www . popia . ft . uam. es / AMIGA 
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2.3 HBT (Han) 

HBT (Hierarchical Bound-Tracing algorithm; Han et al. 2012) is 
a tracking (sub)halo finder. Isolated haloes are first identified with 
a standard Friends-of-Friends (FoF) algorithm (Davis et al. 1985). 
Within each isolated halo, the self-bound part is defined as a cen¬ 
tral subhalo. Starting from the highest redshift, subhaloes are then 
tracked down to later snapshots to link to their descendent haloes, 
by finding host haloes for the progenitor particles. When two or 
more subhaloes are linked to a common descendent halo, we com¬ 
pare the current self-bound mass of the progenitor subhaloes, and 
define their self-bound remnants, except the most massive remnant, 
as satellite subhaloes. The current central subhalo is re-defined to 
be the self-bound part out of all the particles in the host exclud¬ 
ing satellite particles, while its progenitor is defined as the one that 
produced the most massive remnant. The tracking process is then 
continued for all the subhaloes including central and satellites down 
to the final output of the simulation. The position and velocity for 
HBT subhaloes are defined using the 25% of particles with the low¬ 
est local potential energy. 


2.4 Rockstar (Behroozi & Mao) 


The Rockstar halo finder- (Behroozi, Wechsler & Wu 2013) 
adaptively shrinks phase-space isodensity contours to identify 
peaks in phase-space density. Particles within an isodensity con¬ 
tour that contains only one peak are grouped into a single halo 
(or subhalo); the halo’s position and velocity are average values 
for particles near the phase-space peak (typically, within 0. IRyii)- 
When an isodensity contour contains multiple peaks, particles are 
assigned to the closest halo in phase space, determined by the met¬ 
ric d{h,p)\ 
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where h is the halo, p is the particle, CT,, is the halo’s current velocity 
dispersion, Vtnax is its current maximum circular velocity, and p,.,> is 
the virial overdensity from Bryan & Norman (1998). Because par¬ 
ticles are assigned to haloes before the final masses of the haloes 
are known, using Vmax and Ov (which are both consistently mea¬ 
sured even deep inside the halo potential well) improves particle 
assignment stability. When two haloes overlap, the halo with the 
larger number of assigned particles is generally assumed to be the 
host halo. However, in cases where two haloes are within a factor 
of 0.6 in Vmax. information on which halo was the host halo at the 
previous timestep is used to determine which halo will be labelled 
the host halo at the current timestep. 


2.5 SubFind (Muldrew & Srisawat) 

SubFind (Springel et al. 2001) identifies gravitationally bound, 
locally overdense regions in a halo. Initially, a FoF finder with link¬ 
ing length b is used to identify haloes to be processed by SubFind. 
The density of the particles within these haloes is then estimated in 
an SPH-like (Smoothed Particle Hydrodynamics) fashion using an 
adaptive kernel interpolation with Ardens neighbours within the full 
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volume. Locally overdense regions are identified by considering 
each particle in order of density and searching for saddle points us¬ 
ing the Nngb nearest neighbours. Particles with a higher density than 
their neighbours are used to define new candidate subhaloes. Par¬ 
ticles with neighbours that are of higher density, and are attached 
to a single substructure, become members of that substructure. Fi¬ 
nally, particles with denser neighbours that are attached to two dif¬ 
ferent substructures are considered saddle points. These candidate 
subhaloes are then iteratively tested for self-boundedness. Start¬ 
ing with the lowest-density saddle point, a hierarchy of substruc¬ 
ture is determined. Subhaloes are defined as self-bound structures 
enclosed within an isodensity contour passing through the saddle 
point and containing a minimum of Nngb particles. Particles that are 
not assigned to any substructure are added to the ‘background halo.’ 
This is the largest subhalo that was found in the FoF halo, which is 
also tested for self-boundedness. For this study we used the param¬ 
eters b = 0.2 and Aijens = A^ngb = 20. A more detailed description 
of SubFind can be found in §4.2 of Springel et al. (2001). 

2.6 VELOCIraptor (Elahi) 

VELOCIraptor^ (Elahi, Thacker & Widrow 2011, a.k.a. 
structure Einder or STE) is a (sub)halo finder that identifies 
objects in a two-step process. Eirst, haloes are identified using a 
FoF algorithm, where candidate haloes identified by a 3DFoF al¬ 
gorithm are pruned of any artificial particle bridges using a 6DFoF 
and the velocity dispersion of the FoF group. The 6DFoF is also 
used to flag major mergers, that is the presence of two (or more) 
large phase-space dense cores in the FoF halo. Here we follow the 
normal convention and treat the smaller object(s) as a subhalo and 
the larger as a host halo. These field objects are then searched for 
substructures by identifying particles that appear to be dynamically 
distinct from the mean halo background, i.e., particles which have 
a local velocity distribution that differs significantly from the 
averaged background halo. These dynamically distinct particles 
are linked with a phase-space FoF algorithm into substructures. 
Since this approach is capable of not only finding subhaloes, but 
the unbound tidal debris surrounding them as well as tidal streams 
from completely disrupted subhaloes, for this analysis we also 
ensure that a group is self-bound. 

In similar-mass mergers, the mean field is an equal combina¬ 
tion of both haloes, thus neither core will contain (many) particles 
that appear locally dynamically distinct. Hence once the cores over¬ 
lap enough in phase-space, the system will no longer be flagged as 
a merger and the smaller core will not necessarily appear as a dy¬ 
namically distinct substructure either. 


3 STATIC MOCK PROFILE TESTS 

Because no common definition exists for the correct prop¬ 
erties of cosmologically-simulated haloes (Knebe et al. 2013b), 
synthetically-generated haloes are one of the few ways to test 
halo finder accuracy. We adopt the spherical mock host halo NFW 
profile described and tested in Muldrew, Pearce & Power (2011), 
which has a mass‘d of M 200 C = 1.04 x IO'^Mq (other parameters 
include R 200 C = 944.0 kpc, = 259.6 kpc, I'^ax = 715 km s^*). 


^ https : //bitbucket . org/pelahi/velociraptor-stf/ 

In this case, the critical density is calculated for a fiat ACDM cosmology 
with Q.M = 0.3, h = 0.73, at z = 0. 
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Figure 1. Tests with overlapping identical mock halo profiles (§3) with a 
1000 km s“* velocity offset, as a function of the distance between the halo 
centres. In these tests, the “host halo” is taken to be the larger of the two 
returned haloes from each halo finder. Top panel: offset between the input 
and recovered position of the host halo profile. Middle panel: recovered 
host halo velocity; note that most of the halo finders here include substruc¬ 
ture when calculating the host velocity. The average velocity of all parti¬ 
cles when the two halo profiles overlap completely is 500 km s“’. Bottom 
panel: Recovered host halo mass compared to the recovered halo mass for 
the largest separation of the two mock haloes (“Isolated halo mass”). 


The halo profile extends to 2.75 Riqqc and is sampled with ~ 1.5 
million particles, each of mass 1.37 x lO^M©. 

We place this halo in the centre of an empty volume and place 
a duplicate copy of the halo with a velocity offset of 1000 km s^* 
and a distance offset between 0 and 2700 kpc. The chosen veloc¬ 
ity offset is typical for cosmological major mergers on first infall 
(Behroozi et al. 2013b); as the haloes’ velocity dispersions are both 
Ov = 754 km s^*, this places their centres at an offset of 1.33<Tv in 
velocity space. Since the profiles are identical, there is no “cor¬ 
rect” host halo or subhalo choice; however, halo finders typically 
tag overlapping/merging haloes as “host” and “subhalo,” respec¬ 
tively. Note that this choice might be arbitrary, but we also decided 
to follow this notion by referring to the larger of the two recovered 
objects as the “host halo” and the smaller as the “subhalo.” Prop¬ 
erties for this “host halo” are therefore calculated on the combined 
profiles, whereas properties of the “subhalo” are (ideally) calcu¬ 
lated on particles from only one of the individual profiles. 

We show results for recovery of the centre, velocity, and mass 
for the host halo in Fig. 1 . The upper-most panel compares the in¬ 
put position to the one returned for the halo tagged as “host” by the 
respective finder. The positional offset is typically always smaller 
than 10 kpc. Flowever, VELOCIraptor calculates halo centres 
as the centre-of-mass of the innermost 10% of particles within the 
halo radius (determined iteratively), which leads to larger offsets 
when the two mock profiles are within 300 kpc. Some differences 
are also notable in the velocity calculations - as presented in the 
middle panel. AHF and SubFind report averaged particle veloc¬ 
ities within the full host halo radius, including particles from the 
subhalo. Flowever, Rockstar uses a velocity measured closer to 
the halo centre (at 0.1 Ryir) and therefore averages fewer particles 
from the subhalo when determining the host halo’s velocity, un¬ 
til the subhalo approaches much closer to the host halo’s centre. 
VELOCIraptor similarly attempts to exclude most substructure 
when calculating host halo velocities. HBT does not include sub¬ 
structure at all when calculating the host halo velocity, and so the 
recovered velocity is relatively independent of the subhalo’s posi¬ 
tion. The same explanation holds for HBT in the bottom panel, 
where the recovered mass is compared against the input host mass. 
We note that for HBT, which requires a sequence of snapshots in 
order to recover subhaloes, the haloes were processed (and particles 
tracked) in order of largest (2700 kpc) to smallest (0 kpc) distance 
separations. 

Fig. 2 shows recovery of the position, velocity, and mass of 
the subhalo, as well as the fraction of contamination from the host 
halo’s particles. As above, comparison is made to the input posi¬ 
tions and velocities of the halo profile closest to the halo finder’s 
recovered subhalo centre. Extremely good agreement (within ~ 1 
kpc) is seen for position recovery, and the velocity recovery is gen¬ 
erally within 10% of the relative host and subhalo velocities. No¬ 
tably, no advantage in accuracy is seen for phase-space algorithms. 
When the haloes are highly overlapping, particle membership is 
determined largely by a velocity cut; yet, since the halo velocity 
distributions also overlap, this results in an asymmetric truncation 
of the subhalo’s velocity distribution. Without a special averaging 
technique, this leads to a systematic positive radial velocity bias. 
The situation is exactly reversed for the position-space algorithms: 
particles within a given radial aperture of the subhalo centre will be 
contaminated with host particles, leading to a systematic negative 
radial velocity bias. That said, this effect is small in magnitude even 
for the worst-case scenario shown here. Naturally, the effect does 
not exist at all for HBT, as particle contamination from the host is 
not an issue. 
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Figure 2. Tests of overlapping identical mock halo profiles (§3) as a function of the distance between the halo centres. In these tests, the “subhalo” is taken 
to be the smaller of the two returned haloes from each halo finder. Top-left panel: errors in recovering the subhalo’s position. Top-right panel: errors in 
recovering the subhalo’s velocity. Bottom-left panel: ratio of the recovered subhalo mass to the input subhalo mass. The amber dashed line shows the mass 
threshold below which it is impossible to recover the true v^ax of the halo (see §3); see, however. Fig. 3. Bottom-right panel: fractional number of particles 
assigned to the subhalo which were originally from the input host halo’s profile. 


In terms of mass recovery, the position-space finders are at 
a severe disadvantage. Ordinarily, position-space finders can col¬ 
lect particles at large radii around substructure density peaks and 
take advantage of the fact that the high-velocity particles belong¬ 
ing to the host halo will be removed in the gravitational unbinding 
stage. In effect, this performs a quasi-phase-space particle selec¬ 
tion. However, since the velocity dispersion of the two test pro¬ 
files is the same, this technique no longer works. The position- 
space finders therefore end up truncating the mass profile of the 
subhalo as soon as Ihe profiles begin lo overlap (Fig. 2). Vmax re¬ 
covery tends to be much better: given an NFW profile which has 
been spherically truncated, v^ax is very insensitive to the amount 
of truncation, as shown in Fig. 3. Yet, Vmax recovery is impossible 
if the halo itself cannot be found, which happens below 300 kpc for 
SubFind and 45 kpc for AHF. These problems for position-space 
finders are also seen, albeit to a lesser extent, for minor mergers 
(Muldrew, Pearce & Power 2011; Knebe et al. 2011). The phase- 
space finders do very well by comparison, with Rockstar and 
VELOCIraptor recovering between 90-95% of the original halo 
mass. When two haloes overlap significantly in phase space, VE¬ 
LOCIRAPTOR treats them as a single halo; however, Rockstar 
continues to recover two haloes as long as the haloes’ innermost 
density peaks are distinguishable. While both VELOCIRAPTOR 


and Rockstar are phase-space algorithms, this difference allows 
Rockstar to recover properties of haloes in major mergers at 
much closer separations. As expected, HBT is able to perfectly 
recover the input halo’s mass. 

Einally, we address the issue of subhalo purity (lower right 
panel of Eig. 2). Rockstar’s subhalo is contaminated at up to the 
10% level by particles from the host halo, which is in agreement 
with the expected fraction of host particles which are L33C7,. off¬ 
set from the host halo’s central velocity. Eor AHF and SubFind, 
the level of conlamination depends on how aggressively Ihey Iry 
to recover the subhalo’s radial profile. SubFind is more conserva¬ 
tive (and purer), with the result that it recovers much less mass as 
compared to AHF. By definition, HBT does not have any purity 
issues; particles initially assigned to the host are never allowed to 
be assigned to the subhalo. 


4 DYNAMIC MOCK INFALL TESTS 

We next consider a more realistic, dynamically simulated test. The 
initial conditions are two identical mock haloes (the same as de¬ 
scribed in §3) placed at an initial separation of 2700 kpc. We con¬ 
sider two initial velocities for the haloes. In the “Freefall” test, the 
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Figure 3. A comparison of fractional mass loss to fractional v^ax loss in the 
recovered mock subhalo properties, compared to the theoretical expectation 
from the NFW profile used. 


haloes are released from rest. In the more cosmologically realistic 
“Dynamic” test, the first halo begins at rest, and the second halo is 
given a 1000 km s^* velocity offset, aimed toward a point offset 
140 kpc from the centre of the first halo. These initial conditions 
(taken to occur at a = 1) were simulated forward in time to a = 1.2 
by PKDGRAV2 (Stadel et al. 2009), including background cos¬ 
mological expansion according to a flat, ACDM cosmology with 
SIm = 0.3 and h = 0.73 at z = 0; the assumed force resolution was 
6.8 kpc. Because an analytic solution for the merger of two haloes 
is not known, the main purpose of this test was to check the consis¬ 
tency over time of the returned halo properties. 

The positions and radii (with masses being proportional to 
the third power of the radii) of haloes in both tests are shown as 
a time-series plot in Fig. 4. As in the static tests, the position- 
space halo finders show marked artificial mass loss as the haloes ap¬ 
proach each other. In addition, these finders show a dramatic “flip¬ 
flopping” between which halo is assigned to be the host halo and 
which is assigned to be the subhalo. To be fair, the extreme sym¬ 
metry of the tests means that all halo finders which do not include 
some temporal information in deciding host and subhalo relation¬ 
ships will show similar behaviour. However, this behaviour is also 
seen in cosmological simulations (§5; see also Tweed et al. 2009; 
Han et al. 2012; Behroozi, Wechsler & Wu 2013; Behroozi et al. 
2013b; Srisawatet al. 2013). Since Rockstar includes temporal 
information for host-subhalo assignments, it does not show flip¬ 
flopping; similarly, HBT is immune. 

One other feature is evident in the “Freefall” test. In ROCK- 
STAR’s results, the subhalo appears to first grow (a < 1.155) and 
then to lose mass again (a > 1.155). This arises because the cores 
of the merging haloes reach the centre of the potential well and 
begin to orbit rapidly, whereas the remainder of the material has a 
much slower orbital period. When the subhalo core is in-phase with 
the velocity of the outer remnants, the recovered mass increases; 
when the subhalo core is out-of-phase, the remnant is assigned to 
the host halo (Behroozi, Wechsler & Wu 2013). This results in a 
periodic oscillation of the recovered halo masses, which is a fun¬ 
damental limit of the phase-space algorithm employed. This prob¬ 
lem is avoided in VELOClRAPXOR because it no longer separates 
haloes once they begin to overlap significantly in phase space. 


Time Since Simulation Start [Gyr] 
0.8 0.9 1.0 1.1 



a 


Time Since Simulation Start [Gyr] 
0.8 0.9 1.0 1.1 



a 


Figure 4. Time series plots for dynamically simulated twin mock haloes; 
the haloes were released from rest at a separation of 2700 kpc in the top 
panel, and given a more realistic initial velocity offset of 1000 km s“’ in the 
bottom panel (see §4). Each circle corresponds to a halo at the given snap¬ 
shot; circle radii correspond to the halo radius—i.e., they are proportional 
to the cube root of the halo mass—and the bold circle corresponds to the 
host halo. The snapshots are equally spaced in scale factor from a = 1.132 
to a = 1.165, about 420 Myr, and cover a close interaction between the 
two haloes (a = 1.15 and a = 1.16 for the freefall and velocity offset cases, 
respectively). Each tick mark on the bottom axis corresponds to a separate 
snapshot. Results from each halo finder have been spatially offset for clarity. 


5 COSMOLOGICAL TESTS 
5.1 Simulation 

We make use of a simulation described in Srisawat et al. (2013) and 
used for several studies emerging out of the Sussing Merger Trees 
comparison project. This simulation used the GADGET-3 code 
(an improved version of the code presented in Springel 2005) to 
simulate 270^ particles (mass resolution: 1.32 x IO^Mq) in aperi¬ 
odic box with side length 88.8 Mpc. The adopted initial conditions 
were taken from the WMAP-7+BAO+Ho best-fit (Komatsu et al. 
2011); i.e., a flat, ACDM cosmology with parameters Q.,„ = 0.272, 
O-h = 0.0455, h = 0.704, = 0.967, and CTs = 0.810. All the halo 

finders in this project analyzed this simulation as part of Avila et al. 
(2014). 
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Figure 5. A snapshot of the cosmological simulation described in §5, show¬ 
ing the major merger we have selected for analysis (red ellipse). The spheres 
in this image correspond to the locations and radii (Riooc) of haloes returned 
by the Rockstar halo finder. Full movies of the merger process for all halo 
finders are available online.^ 


5.2 Individual Test Case 

We first examine a test case with an individual major merger, se¬ 
lected from the history of the second largest halo in the box at 
z = 0.^ A snapshot of the merger (mass ratio = 1:1.8) is shown 
in Fig. 5, and movies of the halo catalogues returned by all finders 
are available online.® The returned haloes, pruned to exclude all but 
the host halo and the merging halo, are shown in the top panel of 
Fig. 6. The results of §3 and §4 are instructive, as many of the same 
findings apply. As in those sections, the position-space halo finders 
struggle to recover the masses of the subhalo, and show clear ex¬ 
changes of host and subhalo relationships (e.g., AHF at a = 0.78, 
and SubFind at a = 0.76). Additionally, Rockstar shows sub¬ 
stantial variation in the subhalo mass when it passes close to the 
centre; due to the coarser time resolution, this effect is dramatically 
exaggerated in comparison to the dynamic merger test in Fig. 4, but 
the explanation is the same (see §4). Tree building algorithms can 
detect and repair this variability to some extent (e.g., Behroozi et al. 
2013b), regardless of the halo finder used. 

HBT, which had shown exemplary performance for the static 
and dynamic tests (§3 and §4), nonetheless has some issues with 
mass recovery in the cosmological test. In Fig. 6, it is clear at the 
beginning (a = 0.64) that all halo finders except for HBT find 
that the lower halo (corresponding to the leftmost halo in Fig. 
5) has a larger mass than the upper halo. At this snapshot, the 
radii of the two merging haloes are barely touching—so it is not 
an issue of distinguishing between host and subhalo particle as¬ 
signment. Instead, haloes which are major mergers tend to accrete 
mass up to (and sometimes within) the virial radius of the larger 
halo (Behroozi et al. 2014). Because the corresponding friends-of- 
friends groups “bridge” well before the haloes come into contact 
(see, e.g., Klypin, Trujillo-Gomez & Primack 2011), HBT limits 


^ The largest halo is still undergoing a major merger at z = 0, so it is not 
possible to have a clean “before” and “after” comparison. 

® http : //slao . Stanford . edu/-behroozi/MM_Movies/ 
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Figure 6. Merger of the two largest haloes shown in Fig. 5. Top panel: 
time series plots of the merger, as returned by the different halo finders (as 
in Fig. 4). Each circle corresponds to a halo at the given snapshot; circle 
radii ar e proportional to the cube root of the halo mass, and the bold circle 
corresponds to the host halo. Results from each halo finder have been spa¬ 
tially offset for clarity. Bottom panel: mass growth history of the eventual 
host halo, normalized to the host halo mass at the last snapshot. The snap¬ 
shots are equally spaced in scale factor from a = 0.64 to a = 0.87, covering 
a range of 3.69 Gyr. Each tick mark on the bottom axis corresponds to a 
separate snapshot. 


the growth of the lower halo by allowing it only to consume parti¬ 
cles belonging to its original FoF group. 

HBT also highlights an important issue with the definition of 
a halo. Unlike the other halo finders, HBT continues to track the 
smaller halo as a subhalo until z = 0. However, at this point, HBT 
finds that the subhalo’s position is within 9 kpc of the host, on the 
same order as the force resolution of the simulation (7 kpc). In ad¬ 
dition, the velocity-space offset is only 100 km s^^ from the host, 
which is very small compared to the velocity dispersion of the host 
(900 km s^*). Since the mass of HBT’s subhalo at z = 0 is 20% 
of the host mass, no algorithm can robustly distinguish the subhalo 
particles from the host halo particles using phase space informa¬ 
tion alone. While there is no doubt that HBT’s subhalo at z = 0 
is a self-bound structure, it is not clear that all applications would 
wish to treat it as separate from its host (see also discussion in §6), 
especially semi-analytical/abundance matching models and merger 
rate calculations. 

In the lower panel of Fig. 6, we show the mass growth of the 
host halo. In mergers, the sum of the two progenitor haloes’ masses 
can easily exceed the final halo mass after the merger; this is es- 
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pecially true in major mergers. Major mergers dramatically raise 
the velocity dispersion in the resulting merged halo, with the re¬ 
sult that many particles are unbound and that many of the remain¬ 
der spend extended amounts of time orbiting beyond the halo ra¬ 
dius (Anderhalden & Diemand 2011; Behroozi, Loeb & Wechsler 
2013). As a consequence, the mass of the host increases by the mass 
of the smaller halo (a = 0.71), but then rapidly falls as the merger 
raises the velocity dispersion. At the end {a = 0.87), the host mass 
only increases by 15% compared to its original value. HBT again 
is the exception here as the two haloes remain clearly separated 
until a = 0.80, which is when the two centres remain at the same 
position (Fig. 6, upper panel). The subhalo then starts to transfer its 
mass to the host, and the end-mass of the merger approaches the 
result of the other finders. 


5.3 Incidence of Major Mergers 

As a gauge of the importance of addressing the issues discussed in 
this paper. Fig. 7 shows the incidence of major mergers as a func¬ 
tion of halo mass and redshift.^ We define a “major merger” in this 
section to be a subhalo with v^ax at least 70% of that of its host. 
Because Vmax scales as the cube root of mass, this Vmax threshold 
corresponds approximately to a mass ratio threshold of 1:2.9. If 
we had instead defined “major merger” in terms of a subhalo mass 
ratio of 1:3, the position-space halo finders would give very differ¬ 
ent results from the others (Fig. 7, right-hand panels) because they 
often cannot recover full mass profiles for massive subhaloes (§3). 
Indeed, due to SubFind’s conservative approach for assigning par¬ 
ticle membership, it finds a factor of over 30 fewer major mergers 
when using the mass definition as compared to the v^ax definition. 
As shown in Fig. 7, the phase-space and temporal halo finders yield 
similar incidences for the two definitions. 

At z = 0, with the Vmax definition, there is modest agreement 
between AHF, Rockstar, VELOCIraptor, and SubFind that 
between 2-6% of haloes are experiencing a major merger; larger 
haloes are slightly more likely to be undergoing a major merger, 
as they have later formation times (Wechsler et al. 2002). As noted 
in §5.2, HBT tracks major mergers for significantly longer than 
the other finders, with the result that its major merger incidence is 
elevated with respect to the others. 

At z = 2, there is significant (1 dex) disagreement in the in¬ 
cidence of major mergers, as well as in the change in incidence 
compared to z = 0. For IO'^Mq host haloes 1000 particles), 
SubFind finds lower incidences of major mergers at z = 2 than at 
z = 0 (by 0.4 dex), whereas AHF and HBT find similar incidences 
at z = 2 as at z = 0, and Rockstar and VELOCIRAPTOR find 
higher incidences at z = 2 than at z = 0 (by 0.4 dex and 0.8 dex, 
respectively). Given the variance even between similar algorithm 
classes, this likely reflects significant disagreements over where to 
truncate subhalo mass profiles (as in the static tests in §3). We note 
that there is apparent convergence in Eig. 7 for the high-mass major 
merger fraction, but this is illusory: there are not enough haloes in 
the simulation volume to determine whether or not the halo find¬ 
ers agree above a mass of fixlO^^M©. AHF shows a significant 
rise in the merger fraction as a function of mass at z = 2, which 
may suggest that its algorithm is sensitive to resolution effects. All 


^ Note that this is separate from the frequency of major mergers, discussed 
in the introduction. The incidence is the frequency times the average length 
of time that the merging subhalo remains distinct. 


other finders follow the theoretical expectation of more self-similar 
incidences as a function of mass. 

We have also briefly investigated the presence of host-subhalo 
swaps in the merger histories of haloes at z = 0, using the merger 
trees generated for Avila et al. (2014). We find incidences of 0- 
20%, which have some dependence on halo mass and halo finder. 
However, we find that the merger tree algorithm has a much 
stronger influence on the number of host-subhalo swaps than any 
other variable (see also Srisawat et al. 2013). Merger tree algo¬ 
rithms that use all particles contained within haloes (i.e., includ¬ 
ing substructure) to match haloes across timesteps are practically 
immune to host-subhalo swaps (e.g., MergerTree applied to 
AHF; Knollmann & Knebe 2009), at the expense of discontinuities 
in halo positions and velocities (Srisawat et al. 2013). The same 
applies for algorithms which explicitly try to match halo proper¬ 
ties across timesteps (e.g., CONSISTENT TREES; Behroozi et al. 
2013b). Algorithms that use only uniquely-assigned particles (i.e., 
excluding substructure) are much more vulnerable to host-subhalo 
swaps, although they do maintain more consistent halo positions 
and velocities. Due to the strong dependence on the merger tree al¬ 
gorithm, we postpone further analysis of this to a future paper on 
the effect of merger tree algorithms on major mergers. 


6 DISCUSSION 

As we have shown in isolated mergers (§3 and §4), position-space 
finders cannot accurately recover subhalo masses in major mergers, 
although recovery of Vmax is not as severely affected. As shown in 
§4, phase-space finders may perform better when the merger is in 
its early stages, but also have problems when the merger is near¬ 
ing completion. Finally, while temporal algorithms can track halo 
masses well in the final stages of the merger, the initial masses of 
the merging haloes may not agree with other approaches (§5.2). It 
would therefore seem that some kind of hybrid approach is neces¬ 
sary for best accuracy—e.g., phase space when the merging halo 
can still reasonably gain mass, and temporal tracking once the 
merging halo is deep within its host. Alternately, information from 
temporal tracking and phase-space could be weighted in a com¬ 
bined metric for particle assignment. 

As discussed in §5.2, temporal tracking of particles raises the 
issue of when two haloes should be considered to have merged. 
For smaller mergers, this question is less ambiguous, because tidal 
stripping removes much of the mass before the merging halo core 
can sink to the centre. However, in major mergers, the dynamical 
friction is such that the merging halo reaches the centre of the host 
halo before tidal stripping removes all the merging halo’s mass. 
When the centres of the two haloes meet in position and veloc¬ 
ity space, tidal “stripping” ceases to be well defined; instead, the 
profiles of the merging halo and the host halo gradually align un¬ 
til they are indistinguishable. While some applications may benefit 
from continuing to track the merging subhalo past this point (e.g., 
recovery of tidal streams; Elahi et al. 2013), many others may wish 
for such subhaloes to be removed or flagged. A reasonable criterion 
for removal may be when the phase-space ellipse containing the in¬ 
nermost N particles of the subhalo (where N ~ 30) also contains 
more than N particles from the host or other subhaloes. 

Section §5.3 shows that the fraction of haloes undergo¬ 
ing major mergers varies by up to a decade across halo find¬ 
ers at z = 2, even when the mergers are tagged by Vmax instead 
of by mass. As the recovery of the haloes’ Vmax prior to the 
merger is very robust, the main interpretation is that the merger 
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Figure 7. Fraction of host haloes undergoing major mergers as a function of mass. In the left panels, a major merger is defined to be where a subhalo has 
at least 70% of the Vn,ax of its host; this enables fairer comparison with the position-space halo finders. In the right panels, a major merger is defined to be 
where a subhalo’s mass is at least 33% that of its host. Regardless of definition, considerable variation in the reported incidence of major mergers exists at 
z = 2. Some halo finders did not recover major mergers in all mass bins, so their curves above are truncated accordingly. The en'or bars show 68% confidence 
intervals for Rockstar only, to avoid excess clutter; they are not shown in bins where only one major merger was found. Host halo masses of 1.3 x IO'^Mq 
(left-hand edge of both panels) correspond to 100-particle haloes. For comparison, the typical collapse mass M* (i.e., where a{M*) = 1.686) is IO'^'^Mq at 
z = 0 and lO^'^M^) at z = 2. 


timescales and mass-loss rates differ enormously between differ¬ 
ent halo finders. Semi-analytical models which depend on major 
mergers for aspects of galaxy formation (e.g., morphology/size 
changes, black hole growth, or star formation triggering) may 
therefore give very different results when using different halo find¬ 
ers. Abundance matching and similar empirical models’ predic¬ 
tions for the clustering of massive galaxies and the stellar mass 
content of clusters (Leauthaud et al. 2012a) may also be affected 
by the choice of halo finder. Until there exists a reliable method 
for determining which subhaloes host visible galaxies (see re¬ 
cent progress in Wetzel & White 2010; Watson, Berlind & Zentner 
2012; Reddick et al. 2013; Watson et al. 2015), this uncertainty will 
persist. 


7 CONCLUSIONS 

We have examined the recovery of host halo and subhalo proper¬ 
ties in major mergers across five different halo finders. Our main 
findings are summarized as: 

(i) Position-space finders recover subhalo positions and veloci¬ 
ties well, as long as they can detect the subhalo (§3). The recovered 


subhalo Vmax is only somewhat biased, but masses are especially 
difficult to recover accurately (§3, 4, 5). 

(ii) Phase-space finders also recover subhalo positions and ve¬ 
locities well, and recover accurate masses in static profile tests (§3). 
However, in dynamic and cosmological tests, phase differences be¬ 
tween the orbiting halo cores and the remaining merger mass can 
cause large, periodic fluctuations in the recovered subhalo masses 
(§4 and 5). This primarily occurs when the merger is nearing com¬ 
pletion. 

(iii) Temporal finders are able to recover positions, velocities, 
and masses extremely well in isolated merger tests (§3 and 4). 
However, the algorithm tested here (HBT) does not always re¬ 
produce the mass growth history of merging haloes prior to the 
haloes merging (§5.2), when position-space and phase-space find¬ 
ers would be expected to give reliable results. In addition, while 
HBT tracks mergers for much longer than other halo finders, the 
resulting subhaloes are not always distinguishable in phase space 
from the host particles (§5.2). 

(iv) The fraction of haloes undergoing major mergers is in rela¬ 
tive agreement across halo finders at z = 0 (except for HBT); how¬ 
ever, there is strong disagreement (> 1 dex) in this fraction by z = 2. 

(v) All halo finders not using some kind of temporal information 
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show host-subhalo relationship swaps (§4 and §5.2). However, the 
merger tree algorithm employed can to a large degree eliminate this 
problem (§5.3). 

These findings suggest caution when interpreting results from the¬ 
oretical models depending on major merger rates at z > 0, and they 
also suggest that future halo finders have ample room for improve¬ 
ment in their treatment of major halo mergers. 
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